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ABSTRACT 


Let  X  =  (X. .)  =  (X'  . . . ,X') ' ,  X!  =  (X...  ,. . .  ,X.  )  be  a  n  *  p  random 
11  _1  ~n  ,1  ll  ip  r 


matrix  with  probability  density  function 


fx(x)  =  |  Z  |  -n/2  q(tr  I_1(x  -en’)  '(x-ep')) 


where  x  e  )(  =  (x:  nxp  matrix  I  rank  of  n  =  p}  ,  y  =  (y..  , . . .  ,y  )  '  e  R  , 

-Ip 

e  =  n-vector  and  I  >  0  (positive  definite).  Assume  that 

q  e  Q  =  (q:  [0,«)  [0,°°),  convex}  and  n  >  p  so  that  X'X  >  0  with 

probability  one.  It  is  proved  that  for  testing  H^:  y  =  0  versus 

2 

the  alternative  H^:  y  4  0,  the  Hotelling’s  T“-test  is  locally  minimax, 
and  for  testing  H^ :  =  0  versus  the  alternative  H1  :  y(,1^  4  0,  the 


1  *  -U) 


appropriate  Hotelling's  T  -test  is  both  UMPI  and  locally  minimax.  In 

the  second  case  y,,..  =  (ylS...,y  )',  P-.  <  P>  and  (v  ,y  ),  Z  are 

-(.I;  I  p^  i  p-j+i  P 

unknown.  The  above  results  generalize  those  of  Giri  and  Kiefer  (Ann. 
Math.  Statist.,  1964)  under  the  assumption  q  'v  normal.  As  a  technical 
tool,  Wijsman’s  representation  theorem  is  used. 


Keywords  and  Phrases:  Elliptically  symmetric  distributions,  Hotelling's 
2 

T  -test,  Hunt-Stein  theorem,  locally  minimax  tests,  maximal  invariant, 
Wijsman's  representation  theorem. 
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0.  Introduction  and  Summar\ 


Let  X=  (X.  .)  =  (X'  . . .  ,X' ) '  ,  X|  =  (X., 
ij  ~1  ~n  xl 


>X_^)  be  a  n*p  random  matrix 


with  probability  density  function 


fY(x)  =  |Efn/2q(tr  r1(x-eu,),(x-ey'))  (1) 

where  x  e  X  =  {x:  nxp  matrix  |  rank  of  x  =  p}  ,  y  =  (y^, . .  .  ,y^)  '  s  R^, 
e  ■  (1, . .  .  ,1)  '  n-vector  and  I>0  (positive  definite  pxp).  We  shall  assume  that 
q£Q  =  {q:  [0,°°)  [0  ,°°) ,  convex}  •  This  is  a  subclass  of  probability  density 

functions  which  are  left  0(n)  (nx  n  orthogonal  matrices)  orthogonally  invariant 
distributions  about  ey'  and  is  also  a  subclass  of  elliptically  symmetric 
distributions  about  ey'  with  scale  matrix  E.  We  shall  assume  throughout  that 
n  >  p  so  that  X'X>0  with  probability  one  (see  Giri  (1977)). 

We  shall  write  for  any  p-vector  b=(b,,...,b  )'  =  (&»  ,b')  with 

1  P  ~(1)  ~(2) 

(l)*^l****’bp  ^ '  ’  ^(2)==^*3p+l’’’’’^p^'  an<^  S[i]  =  (b^,...,b^)'  and  for  any 


2 


px  p  matrix  A=  (a.  .)  = 


A(ll)  A ( 12 ) 
A ( 21 )  A(22)J 


with  A(11)  »A(22) »  P^Pj.  «nd  P2*P2  sub¬ 


matrices  respectively,  satisfying  p1+p2  =  p.  Also  we  shall  write 


A[i]  = 


11' 


*’ali 


;aii’--*’aii 


the  ix i  left-hand  corner  submatrix  of  A.  Denote  by 


1  n  n 

X  =  -  Tx.  ,  S  =  J(X.-X)(X.-X)'. 

~  n  ,  ~i  ~  -i  ~ 


1  1 

The  assumption  of  multivariate  normality  i.e.  when 


q (tr  E  1(x  -  ey* ) ' (x  -  ey' ) ) 

-np 

=  (2  IT)  2  exp{-  tr  E_1(x  -  ey' ) '  (x  -  euf)  } 


(2) 


leads  us  to  derive  optimum  test  procedures  for  testing  problems  concerning 
U  and  E.  However  it  has  been  established  that  (see  for  example  Kariya  and 
Sinha  (1984)  and  the  references  contained  therein)  optimum  procedures  can 
also  be  derived  by  replacing  the  multinormality  assumption  by  one  closely  re¬ 
lated  to  it,  namely,  the  class  of  spherically  symmetric  distributions  in  the 
case  of  a  random  vector  X  (pxl)  and  the  class  of  elliptically  symmetric  dis¬ 
tributions  in  the  case  of  an  nxp  random  matrix  X.  These  optimum  procedures  in 
turn  involve  the  robustness  of  the  distribution  of  the  test  statistic  under 
the  null  hypothesis  (null  robustness),  the  robustness  of  the  distribution  of 
the  test  statistic  under  the  alternatives  (non-null  robustness)  and  the  ro¬ 
bustness  of  the  optimum  properties  of  the  test  procedures  (optimality  robust¬ 
ness).  We  shall  call  a  test  robust  in  any  one  sense  if  an  optimality  property 
which  the  test  enjoys  can  be  extended  to  a  class  of  distributions  including 
the  distribution  under  which  the  optimality  holds. 


3 


(X,S;  y,Z)  -  (gX,  gSg';  gy,  gig'),  geG£(p). 

2  _i_ 

A  maximal  invariant  in  the  space  (X,S)  under  Gf (p)  is  T  =nX'S  X  or  equivalently 

^  _ 2  2 

R“  =  nX'  (S  +  nXX' )  X  =  T  /(1+T  ).  A  corresponding  maximal  invariant  in  the 

2  -1 

parametric  space  of  (y,I)  is  5  =  ny '  I  ’ i .  Kariya  (1981)  has  proved  that  the  Hotelling' 

2  2  2 
T  test  which  rejects  whenever  T  >C  or  equivalently  R  _>C,  where  C  is  a 

constant  depending  on  the  level  a  of  the  test,  is  uniformly  most  powerful  in- 

2 

variant  (.UMPI),  whatever  q  e  Q.  The  distribution  of  T  under  is  the  same 
2 

as  that  of  T  under  the  multivariate  normal  set-up  (Kariya,  1982).  Needless 

to  mention  that  the  multivariate  normal  distribution  belongs  to  the  family 

2 

given  in  (1).  We  shall  show  in  section  1  that  the  Hotelling's  T  test  is  lo¬ 
cally  minimax  in  the  sense  of  Giri  and  Kiefer  (1964)  for  testing  H^^  against 
H^  in  (1)  as  '“-O.  In  the  multivariate  normal  setup,  this  result  is  proved 
:n  Giri  and  Kiefer  (1934). 

Let  be  the  group  of  translations  such  that  t^  ?  translates  the  last 
p0  components  of  each  X  ,  i=l,...,n  and  let  G  be  the  multiplicative  group  of 
p-p  nonsingular  matrices  of  the  form 


4 


where  is  the  upper  left  hand  corner  p^xp^  submatrix  of  g.  The  second 

problem  remains  invariant  under  the  affine  group  (G,T^)  such  that  for  g  G, 


t  -  T 
1  '  1 


( 8 ,  t ^ ) X  ~  b i ,  i  1, . . .  ,n. 

A  maximal  invariant  under  (G,T^)  is  given  by  (see  Giri  (1977)) 

*1  =  n?(l)(S(ll)  +r?(l)?(l))  ~(1) 

=  n?(l)S(ll)?(ll)/(1 +n?(l)S(ll)?(l))  (5) 

A  corresponding  maximal  invariant  in  the  parametric  space,  under  the  induced 
group  of  transformations,  is 

^1  =  n^Cl) E (11)~ (1) *  (6) 

For  any  invariant  test  under  (G,T^),  the  second  problem  reduces  to  testing 

-  2  -  2 
^20:  ^1  against  the  alternatives  >  s'ria^l  show  in  section 

-2 

2  that  the  test  which  rejects  whenever  _>C  is  UMPI  and  locally  minimax 

_2 

as  5^-*0  for  (1).  In  the  multivariate  normal  setup  this  test  has  been  proved 

_2 

to  be  UMPI  and  locally  minimax  as  (Giri  and  Kiefer  (1964)). 


1.  Locally  minimax  test  for  problem  1 

The  theory  of  locally  minimax  tests  has  been  developed  in  Giri  and  Kiefer 

.  2 

(1964).  We  refer  to  this  paper  for  details.  For  each  (h  , rj )  in  the  parametric 


space  2,  where  a  >0  and  n  is  of  arbitrary  dimension  and  its  range  may  depend 
2  -2 

on  6  ,  let  p(x;  ,n)  be  a  probability  density  function  on  (X,A)  with  respect 


5 


to  some  a-finite  measure.  Suppose  that  we  are  interested  in  testing  at  level 

2  *  2 
(0  <  o  <  1)  the  hypothesis  H^:  6  =0  against  the  alternatives  H^:  5  =  A 

where  A  is  a  positive  constant.  For  fixed  a,  consider  the  critical  region  of 

the  form 


where  0  c^.  <r(-\)<c^<M  for  A  sufficiently  small  and  g(A)=o(l)  and  B(x,A)= 

(h(A))  uniformly  in  x.  If  U  satisfies  (8)  and  (9)  and  if  for  sufficiently 

small  A  there  exist  and  £  ^  satisfying  (10)  then  is  locally  minimax  for 
2  2 

testing  H^:  >  =0  against  the  alternatives  HJ^:  $  =  A  (specified)  as  A  -*•  0 . 

It  is  wellknown  that  (see  for  example  Giri,  Kiefer  and  Stein  (1963), 

Giri,  Kiefer  (1964))  the  Hunt-Stein  Theorem  cannot  be  applied  for  the  group 
G?(p)  with  p  2.  However  this  theorem  does  apply  for  the  smaller  group 


6 


G  =  (pxp  nonsingular  lower  triangular  matrices} 

2 

with  p >_2.  Thus  for  each  5  ,  there  is  a  level  a  test  which  is  invariant  under 

G^  (see  Lehmann  (1959),  p.  225)  and  which  minimizes  among  all  level  o.  tests, 

2 

the  minimum  power  under  In  the  place  of  R  under  G^(p),  we  obtain  a 

2  2 

p-dimensional  vector  R=  (R_,...,R  ) '  as  a  maximal  invariant  under  G^  and 

-  1  p  T 

R  is  defined  by 

j1Rj-n?ii]<s[iii+n  ?[i]?ii])_1B[i] 

i  -  1 , • • •  ,P 

2  P  2  2 

with  R^j>0  and  )  R.  =  R  .  A  corresponding  maximal  invariant  in  the  parametric 

„2  9  2  2 

space  under  the  induced  group  is  0  =  (0.  , . . .  ,5“) '  and  are  given  bv 

—  Ip  In 


r  ?2  _  ,  _-l 

J  ~  [ 1 ]  [ XI ] ~ [ X  ]  * 


,2  2  .2  .2 

with  ;  i  = 5  .  The  nuisance  parameter  in  this  reduced  setup  is 

1  c2  -2 

5  o 

n  -  (-r  -£)’  =  (n, , . . .  ,n  )'.  (13) 

52  62  1  P 

Since  the  distribution  of  R  depends  on  Q  only  through  5  we  may  put  Z =  I  and 
redefine  *  n  L  ~  ±  =  (b^, .  • .  ,5  ) ' .  Using  Stein's  theorem  (1956)  or  Wijsman's 
representation  theorem  (1967),  the  ratio  of  the  probability  densities  of  R, 
under  H'  and  under  H  ,  is  given  by  (for  g  G  ) 


2  ,  2 


rat io= 


Pi (gx,gsg' )H(g“ .)  dg 

»  ^  i  -1-  A- 


j  P0(gx,gsg' )H(g"  )  2  dg 
JG^  u  ~  1  11 


where 


,i., .i-li  ^hih  .^  p  i ^ up ^"i.1  r  rrr»rv»:* r* rrr ».  i. 

7 

pl(x,s) 

PQ(x,s) 


=  |  s  |  2  q  ( t  r  { ( s+nxx f )  -  2nx '  u4-S  } ) 


=  [ s | “  q (tr ( (s+nxx’ ) ) . 


Lee  A  be  a  matrix  belonging  to  G  such  that  A(s+nxx')A'  =  Ip-  Then 
A'A=  (s+nxx')  ^  =  s  1  -  ns  ^xx's  3(l+nx's  \x)  \  so  that  nx'A'Ax  = 


nx's  ^x(l  +  nx's  3x)  3  =  ^R7  =  R2. 


Since  A r  .  .  .  (s  r  .  .  ,  +  nx  ,  .  ,  x!  )a\  .  .  ,  =  Ij  » 
[ll]  [n]  ~[i]~[i]  [n]  1 


we 


obtain  nx! . ,A! . . ,Ar . . ,x, . ,  =  so  that 

~  L i  ]  In]  [ ii ] ~  [ i ]  r  j 


n  Ax  =  (R  , . . . ,R  ) '  =2 
1  P 


(15) 


where  v  !  .  n  v  , .  ,  -  / i\ . 

~  [  i  ]  ~  [  i  ]  “  j 


r  2  -1 

■R..  Writing  gA  =  g  we  get  from  (14) 


n-i 


ratio  = 


q(tr(gg'  -  2  a'gy  +  52))  IT  (g2  )  2  dg 
gt  -  i=l  11 


r 


n-i 


q ( tr (gg ' )  H(g2  )  2  dg 
CT  1 


(16) 


where  i  =  (.5,  '  ) '  . 

1  P 

Let  us  assume  that  q  is  thrice  continuously  differentiable.  Writing 

(i),  ^  dj.  ,  „  . 

q  (x)  =  — ^  we  obtain 

dx1 


q  ( t  r  ( gg '  -  2  :i'g£  +  5  )) 


=  q(tr  gg')  +q('1)  (tr  gg’ )  (-2tr  (a*  gy)-K\) 

(2)  , 

+  —  (tr  gg')(-2  tr  ff'gxHA)" 


(17) 


+  q(3)(Z)(-2  tr  {J  gx)+A)3 


where  Z  =  tr  gg '  +  (1-a)  (-2  tr  (u'  g^)+A)  ,  0<_^_<1.  Let 

n-i 

P  2  ? 

D  =  q(tr  gg')  IT(g  )  “  dg. 

jgt  1  11 


(18) 


n-i 


(i)  p  2  2 

Since  the  measures  q  (tr  gg')  IKg^)  dg,  i  =  1,2,3  are  invariant  under 


the  change  of  sign  of  g  to  -g,  we  get 


n-x 


J 


Ci')  p  2  2 

tr(a'g£)qk  '(tr  gg')  II  (g  .)  dg  =  0 

gt  i-i 


and 


n-i 


f  C  \  P  2  2 

J  g^goiV  ;(tr  gg’>  n  (g . , )  dg  =  o 


3ij  6  2k 


li 


i=l 

if  i  ^  i,  j  ^  k. 


Thus  from  (16 )  —  ( 2 0 )  we  get 


n-i 


if  / 1  \  P  o  2 

ratio  =  l+£  I  qV  '  (tr  gg^IUgf  )  dg 

jgt  1 


n-i 


+  f  |  (tr  a'g^)2q^2)(tr  gg^JUg^)  2  dg 


n-i 


2D 


q^2) (tr  gg')n(g2  )  2  dg 
GI  1 


+  ^  j  (-2  tra'gy  +  A)3  q (3)  (z)II(g2i)  2  dg. 


n-i 


(19) 


(20) 


(21) 


The  first  integral  in  (21)  is  a  finite  constant  8^.  To  evaluate  the  second 


integral  in  (21)  we  first  note  that 


tr  a’g>  =  l  r  [  l  5  g  +  6  g  ], 

j  J  i>j  J  J  JJ 


(22) 


From  (20)  and  (22)  the  second  integral  can  be  v ritten  as 


n-i 


•r  2  v  ^2  2  ■  x2  2  i  (2).  ,  P,  2.2 

(Z  r,  )  (5igii  +  5i8ii)9  (tr  gg')IUg„)  dg. 

-t  J  ■[>  j  J  JJJ 


G_  i 
T 


(23) 
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Let  us  nov:  evaluate  the  integral 


Define 


f  2  (2 'l  p  ?  2 

i  =  |  g  ,qU)(u  gg’)n(g;  )  dg. 
jgt  11  1  11 


L  =  tr(gg’) , 

Si  =  gii/L  5  i=1>2 . P  5 


Si+l,i/L’  1  1’*** ’P_1  ; 


6  i , •  2 . , «  • f  1 

p+p-l+i  °i+2#i’  r 


■  gpi/L- 


Write 


K  -  L  q  '  (tr  gg')dg 


r  i 

N  =  I  L 


(— ) 
v  2  ; 


q(2)(L)dg. 


(2)  f 

Since  ^ -  is  a  spherical  density  of  g_’s,  L  and  e=  (e^, . . .  ,e^^+^ ) ' 

11  2 

are  independent  and  e  obeys  a  Dirichlet  distribution  D (y, . . .  .  From 

Kariya  and  Eaton  (1977)  the  probability  density  function  of  e  is  given  by 


p  (e)  = 


[r(i)]p(p+i)/2 


2  1  i 

n  (e  )^_1  (1  - 

i 


i-1 


Now  using  (24)-(26)  we  get 


n-j  n-i+2 

N  2  2 

I  -  7  E(  n  (e.)  e.  Z  ) 

K  .  ,  .  j  l 

=  N  (n-i+1) 

K  2  ’ 


(27) 


t-”’-'”'- — ■  V  'VV  ". '  V  ■*. 'A  1  ■.  ■  .  *.  ■'  ■„  ^  ■.  1  V1  1  ■*'  «. 
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where 


M  = 


P 

E(n(e.) 

i 


n-i 

2 


). 


Similarly 


n-i 


2  (?)  p  2  ? 

g:4qU;(tr  gg'inCgT.)  2 


ij 


'n 


dg 


P  n-i 

=  E (e  He,  2  ) 
k  1  x 


NM 

=  JE 


wi  th  k 


> 


P  * 


(28) 


It  may  be  verified  that  the  third  and  the  fourth  integrals  in  (21)  are  both 

r2 

of  the  order  o(,«  )  uniformly  in  and  0.  Hence  we  can  write 


ratio  = 


(n-j+l)n2) ) 


+  B(j,n,A) 


(29) 


where  B(^,_n,A)  =  o(A)  uniformly  in  £  and  ru  Thus,  from  (29),  the  equation 

(10)  is  satisfied  by  letting  give  measure  one  to  the  single  point  n_  = 

*  *  * 

while  £•>  give  measure  one  to  the  single  point  0  =  n  =  (n,  )  whose 

IA  —  — *  x  p 

*  « 
j th  coordinate  n.  satisfies 

j 

*  -1  -1  -1 
=  (n-j)  (n-j+1)  p  n(n-p),  j=l,...,p 


so  that 


-  *2  *2  n 

l  n.  +  (n-j+Dn,  =  f.  (30) 

i>j  1  3  P 

2 

The  equation  (8)  follows  from  the  null  robustness  of  *"he  distribution  of  T 
(Kariya  (1982)).  To  show  (9)  we  proceed  as  follows. 


i  o 


For  g€G0,  we  can  write  using  Wijsman's  theorem 


f  2(t2|H'  )  j  q(tr(g§'  “  2a'g^+X)  |gg*  |  2  dg 

f T2(t2|Hio)  r  “ 

q(tr  gg')|gg' |  dg 


To  evaluate  (31)  we  use  the  following  results. 

(i)  Given  x.  =  (x^,...,x  )’  there  exists  an  OeO(p),  the  group  of  p*p 
orthogonal  matrices^  such  that 


Ox  =  (/x' x  ,  0,...,0)'. 


(ii)  We  can  decompose  G^  =  GT*0(p).  Denote 


Hl£ 

•  I  2 


M (dg)  =  | gg ' |  dg,  K  (dh)  =  n(h2  )  2  dh 

1  X 


for  h €  G^  and  x(dO)  the  invariant  measure  on  0(p).  Then  for  g?  G^, 
U(dg)  =  C (dh)  T(dO). 


(iii) 


tr(AOBO')x(dO)  = 


and  [  (tr  OA)  T(dO)  = 

J  0(p) 


r  w 

trA* 

-  r\ 


0  if  k  is  odd 


if  k  =  2. 


Using  the  above  results  we  can  rewrite  (31)  as 


\  r  Cl>  .  P  o  9 
1  +  £  I  qU'(tr  hfi)  n(hf.)  Z  dh 

D  Jgt  1  11 


2  (2) 


2  \  2 


+  ~2~  hUq  (tr  hh')n(h^)  A  dh  + 


P  D  J  G 


=  1  +  \ (K  +  cr  )  +  B(y,n,A) 
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where  B(y_,n,A)  =  o(A)  uniformly  in  £  and  n  and  k,c  are  positive  constants 

2 

(using  (27)  and  (28)).  Hence  for  T  -test,  with  b>G, 

P.  (^)  =  o+  bA  +  g(X,n) 

a  ,n 

with  g(X,H)  =  o(bX). 

Hence  we  get  the  following  theorem. 

2 

Theorem  1.  For  testing  H^  against  H^,  Hotelling's  T  test  is  locally  mini¬ 
max  as  >.->0  for  the  family  of  distributions  given  in  (1). 
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where,  taking  w  =  ^22^1U* 


q(tr  Z(ii)V)  = 


R 


P2n 


q(tr(E^^V  +  ww’))dw. 


The  assumed  convexity  of  q  implies  convexity  of  q. 

To  find  the  ratio  of  the  probability  densities  of  the  maximal  invariant 

_2 

with  respect  to  the  group  (G,T^)  we  can  without  any  loss  of  generality 
consider  only  the  multiplicative  group  G^(p^)which  transforms  each 

Ji(i)*sJi(D'  seGt(pi)- 

Using  Stein's  theorem  or  the  Wijsman's  representation  theorem  we  obtain 

-2 

for  geG^(p^),  the  ratio  of  probability  densities  of  R^  under  and 
H20  aS 


n-p, 


ratio 


Jl,((pr,()88'"2- 81  +6i 

io  =  — - 1 - 


gg' I  2 


dg 


n-pi 


|q(tr  gg') J  gg ' (  2  dg 

JVpx) 


(38) 


2  _2 

where  £  =  (p  r  ..  ,p  ) ' ,  £  =  (r  , . . .  ,r  ) '  satisfying  £'p  =  5  ,  x'x  =  ri  • 
i  P-^  i  p^  -'I  i 

Using  (32)  we  can  write  the  numerator  of  (38)  as 

n-Pl 


q(tr(gg'-2£*'g£*  +  52)|gg'|  2  dg  =  H(rj)  (say) 


with  £*  =  (6^,0, . . . ,0) 1 ,  £*  =  (r^,Q, . . .  ,0) ’ .  Since  |  gg 


n-Pj 

.1  2 


under  the  sign  change,  g-^-g  =  hfG^p^)  we  conclude  that 

n-Pi 


H 


(r  )  =  j q ( tr (gg' -2p*' g(-£*)+?2) j gg ' |  2  dg  =  H(-r  ), 

JG?(Pl) 


dg  is  invariant 


(39) 


Hence  H(r^)  =  xH(r^)  +  (l-t)H(-r^).  Since  q  is  convex  by  assumption,  for 


nj|h-“ 
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<  a  <  1, 


H(r^)  j>  H( (2a-l)r^) .  (40) 

From  C 39 )  —  (40 ) ,  we  conclude  that  H(r^)  is  a  monotonically  increasing  funciton 

-2 

of  r^.  Now  applying  Neyman-Pearson  lemma  we  obtain  the  following  theorem. 

_ 2  2 

Theorem  2.  For  testing  ^q:  6^  =  0  against  the  alternative  "o  >  0  the 

-2 

test  which  rejects  whenever  R  _>  c  is  UHP1  with  respect  to  the  group 
(G,T^)  for  the  family  of  distributions  given  in  (1). 

Moreover,  we  can  prove  in  an  analogous  wayj using  the  results  of  section  1, 
the  following  theorem.  The  details  are  omitted. 

_ 2 

Theorem  3.  For  testing  H^q  against  6^  =  X  (specified)  >  0>  Hotelling' s 

2  -*> 

T  test  which  rejects  H^q  whenever  R~  _>  c  is  locally  minimax  as  X  0. 
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